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Abstract 

An autopilot can be used to provide precise 
control to meet the demanding requirements of 
flight research maneuvers with high-performance 
aircraft. This paper presents the development of 
control laws within the context of flight test 
maneuver requirements. The control laws are 
developed using eigensystem assignment and com- 
mand generator tracking. The eigenvalues and 
eigenvectors are chosen to provide the necessary 
handling qualities, while the command generator 
tracking enables the tracking of a specified state 
during the maneuver. The effectiveness of the 
control laws is illustrated by their application 
to an F-15 aircraft to ensure acceptable aircraft 
performance during a maneuver. 

Introdu ction 

Conventional piloting techniques are often 
inadequate to meet the demanding requirements of 
flight research maneuvers with high-performance 
aircraft. These maneuvers often require precise 
control of onset rates in extreme flight condi- 
tions. Thus the pilot may be trying to control 
an aircraft at high angles of attack and high g's 
while attempting to increase normal acceleration 
at a prescribed rate through a maneuver specified 
at the very limits of accuracy of the cockpit 
i nstruments. 

A new flight test technique 3 was developed at 
the Dryden Flight Research Facility of the NASA 
Ames Research Center (Ames-Dryden) to aid the 
pilot during these maneuvers. The essence of this 
technique is the application of an autopilot to 
provide precise control during the required flight 
test maneuvers. The flight test maneuver auto- 
pilot (FTMAP ) is designed to provide precise, 
repeatable control of a high-performance aircraft 
during certain prescribed maneuvers so that a 
large quantity of data can be obtained in a mini- 
mum of fl ight time. 

The FTMAP can be used for various maneuvers, 
including strai ght-and-level flight, level accel- 
erations and decelerations, pushover pullups, 
excess-thrust windup turns, thrust-limited turns, 
and the "rockinghorse" maneuvers. Each of these 
maneuvers involves tracking certain states of the 
aircraft and holding certain states within pre- 
scribed values, as well as placing constraints on 
the derivatives of the states (for example, the 
excess-thrust windup turn is performed at constant 
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altitude and Mach number with angle of attack 
increasing at a specified rate to achieve the 
final angle of attack). 

This paper discusses the development of the 
FTMAP control laws within the context of flight 
test maneuver requirements. The FTMAP control 
laws are developed using the eigensystem assign- 
ment 2 and command generator tracking. 3 The feed- 
back gain obtained from eigensystem assignment 
ensures the desirable performance of the aircraft. 
The feedforward methodology of command generator 
tracking ensures the tracking of desired states 
during a particular maneuver. The desired closed- 
loop eigenvalues are taken from MIL-F-8785-C,^ and 
each of these eigenvalues is assigned an eigenvec- 
tor that distributes the modal response among the 
state variables and outputs of the system. How- 
ever, each eigenvector is constrained to lie in an 
m-dimensional subspace 5 (m is the number of inde- 
pendent controls), and an element in this subspace 
is selected by finding the best linear projection 
of an unconstrained desired vector in this sub- 
space. A feedback gain matrix, K, is computed 
using the eigenvalues and achievable eigenvectors. 

The feedforward gains are computed by using the 
command generator tracking of Broussard. 3 This 
method would ensure the tracking of a state trajec- 
tory while the aircraft is performing a particular 
maneuver. The tracking objective can be described 
in terms of the controlled variables of the air- 
craft being able to follow the output of a model. 
The input to the model, which is the pilot input, 
is assumed to be constant and the output of the 
model is the desired state trajectory the aircraft 
will follow while undergoing a maneuver. 

The effectiveness of the developed control laws 
is illustrated in this paper by application to an 
F-15 aircraft. 

Eigensystem Synthesis 

Two widely used synthesis techniques of modern 
control theory are the linear quadratic regulator 
design and the modal control theory involving pole 
placement or eigenvalue/eigenvector assignment. 

One of the purposes of feedback control of aircraft 
is to improve or enhance the flying qualities of 
an aircraft. The difficulty in incorporating spec- 
ifications such as damping, natural frequency, 
and decoupling within a quadratic performance 
index makes the eigensystem synthesis procedure a 



promising design alternative. The performance 
specifications can be interpreted in terms of 
desired closed-loop eigenvalues and ei genvectors . 

Moore® and others have shown how feedback can be 
used to place closed-loop eigenvalue and shape 

closed-loop eigenvectors. Cunningham, 7 Andry, 
Shapiro, and Chung, 8 and Sobel and Shapiro 9 have 
successfully demonstrated the use of eigenstruc- 
ture assignment procedure for aircraft control 
system design. 

The handling qualities data base may be used 
to obtain desired pole locations directly. The 
additional design objective of obtaining augmented 
dynamics similar to those obtained in flight leads 
to specifications on eigenvectors or desired mode 
shapes. For example, pitch attitude must be domi- 
nant for the short-period mode, and speed must be 
dominant for the phugoid mode. 

Detailed discussions on eigenspace can be found 
in Ref. 8, but some basic results for controllable 
and observable systems are summarized in the fol- 
lowing discussion. 10 

Consider the system 


x = Ax + Bu 

(1) 

<< 

It 

O 

X 

(2) 


where x e R n , u e R m , and y e R r . A, B, and C are 
matrices of appropriate dimensions. If the system 
is controllable and observable, and the matrices B 
and C are of full rank, the following results 
hold. 


only a few of the components in v^ are actually 

specified. The rest can be arbitrary. To account 
for this, we reorder and partition v- as follows: 

[.,]"< .(^l) (3) 

where 

v* is the specified subvector 

d^ is the vector of unspecified components and 
R 

[...]■> is the reordering operation 
If we let 

,(£)" [( V - Arl ^ Ri Z 1 =ft) z i (4) 

then, as shown in Ref. 5, we can select Z- to best 

approximate v* with v? and by method of orthogonal 
projections obtain 

Z i = (Car 1 L'v* (5) 

As shown by Moore,® the feedback gain K is given by 
K = [Wj w 2 . . . w n ][ Vl v 2 . . . v^' 2 3 4 1 (6) 


1. The positions of maximum (m, r) closed-loop 
eigenvalues can be assigned arbitrarily with the 
stipulation that if is a complex closed-loop 
eigenvalue, its complex conjugate, X*, must also 
be a closed-loop eigenvalue. 

2. The shape of maximum (m, r) eigenvectors can 
be altered. If the shape of a complex eigenvector 
v^ is altered, its complex conjugate, v*, must be 
altered in the same way. 


where w^ is obtained from the relation 

[X i I - A] Vi = Bw. (7) 

Command Generator Tracking 

The development of control laws based on com- 
mand generator tracking is discussed by Broussard. 8 
Briefly, if a model specifying the desired behavior 
of an aircraft is defined by 


3. For each eigenvector whose shape is altered, 
minimum (m, r) eigenvector elements can be chosen 
arbitrari ly. 

4. Attainable eigenvectors must lie in the sub- 
space spanned by the columns of (X..I - A) _1 B of 
dimension m, which is the number of independent 
control variables. A desired eigenvector, v^, 

will, in general, not reside in the prescribed 
subspace and cannot be achieved. The achievable 

eigenvector, v^ , is obtained by orthogonal pro- 
jection of v*j onto the subspace spanned by 
(X i I - A)'*B. It will generally be true that 



( 8 ) 

(9) 


and the variables of the aircraft to be tracked 
are described by 

y t = Hx (10) 


a control law must be found such that y^ - y^. If 

x* and u* are the ideal aircraft state and input, 
respectively, that allow perfect tracking, then 
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x* = S,,x m + 
11 m 

S,,u m 
12 m 

(id 

A m = 0 

u* = S 91 x m + 
21 m 

SrtrtU 
cl m 

(12) 

Bm = I 


The results in Eqs. (11) and (12) assume that the 
input to the model, u m , is a step input. However, 
it does not mean that only a constant command 
input is being tracked. Any command signal that 
can be described as a solution of the differential 
equation forced by a step input (or zero) can be 
used, provided it is augmented to the model state 
and not to the model input. 1 ! The matrices S^, 

S 12’ S 21’ and S 22 are given by 


S,, = SI. ,S. ,A + S1.„C 
11 11 11 m 12 m 


S.„ = S2. ,S, ,B + SJ. „D 
12 11 11 m 12 m 


S 21 = n 21 S ll A m + n 22 C m 


S 22 = S2 21 S ll B m + fl 22 D m 


and 


fl = 


'On 

“l2 


"A 

B~ 

“21 

“22. 

II 

_H 

V 


-1 


(13) 


(14) 


To incorporate the state feedback into the design, 
i f we 1 et 


x = x 
u = u 


x* 

u* 


then 


and 


x = Ax + Bu 


0 = K(x - x*) 


u = u* + u = S„.x m + S„ 0 u m + K(x 
21 m 22 m 


(15) 

(16) 

(17) 

(18) 


x*) 


= S 


21 x m 


+ S„„u + Kx 
22 m 


K(S ll x m + S 12 U m^ 


(19) 


Equation (19) gives the control to be applied to 
the aircraft, depending on the maneuver of the 
aircraft. For example, for constant-level accel- 
eration, the command velocity is to increase 
linearly while the command altitude is held 
constant. The input to the model is constant, 
while the outputs of velocity and altitude are a 
ramp function and zero, respectively. In this 
case. 


Cm = 1 

D m = 0 (20) 

with an appropriate value for u m . In addition, 
Sn = n 12 

s 12 = « 11«12 
S 21 = “22 

s 22 = « 21«12 ( 21 ) 
x * = n 12 x m + nil«12 u m 
u* = « 22 x m + 021«12 u m 


( 22 ) 

(23) 


and 


u = (022 - K0i2)xrn + («21«12 

v 


K0nfti2)u m 


+ Kx 


Feedforward gain 


Example 


(24) 


The feedback and feedforward control technique 
developed by eigensystem and command generator 
tracking synthesis technique is applied for the 
control of an F-15 aircraft. Ames-Dryden's 
detailed nonlinear aerodynamic model of the F-15, 
linearized by trimming the aircraft at the desired 
flight condition and deriving linear models by 
numerical perturbation, is used. The linear model 
is represented in the state equation form as 


where 


x = 


u = 


L t J 


Ax + Bu 


velocity 
angle of attack 
pitch rate 
pitch angle 
altitude 

angle of sideslip 
roll rate 
yaw rate 
roll angle 


aileron deflection 
elevator deflection 
rudder deflection 
throttle 


For the flight condition corresponding to an 
altitude of h = 20,000 ft and a Mach number (M) of 
0.8, the values for A and B matrices of Eq. (1) 
are as follows. 
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--0.0108 

24.1966 

0.0000 



-0.0001 

-1.0942 

1.0000 



-0.0001 

-3.2862 

-2.1922 



0.0000 

0.0000 

1.0000 


A 

= 

0.0000 

-829.5390 

0.0000 



0.0000 

0.0000 

0.0000 



0.0000 

0.0000 

0.0000 



0.0000 

0.0000 

0.0000 



_ 0.0000 

0.0000 

0.0000 



-32.1129 

0.0000 

0.0000 



0.0001 

0.0000 

0.0000 



0.0009 

0.0000 

0.0000 



0.0000 

0.0000 

0.0000 



829.5390 

0.0000 

0.0000 



0.0000 

0.0000 

-0.2337 



0.0000 

0.0000 

-40.0103 



0.0000 

0.0000 

9.0098 



0.0000 

0.0000 

0.0000 



0.0000 

0.0000 

0.0000“ 



0.0000 

0.0000 

0.0000 



0.0000 

0.0000 

0.0000 



0.0000 

0.0000 

0.0000 



0.0000 

0.0000 

0.0000 



0.0358 

-0.9994 

0.0387 



-2.1420 

1.2406 

0.0000 



-0.0340 

-0.6040 

0.0000 



1.0000 

0.0358 

0.0000_ 



- 

0.0000 

-1.0734 

0.0000 0 

.3792“ 



0.0000 

-0.1504 

0.0000 0 

.0000 



0.0000 

-16.1223 ' 

0.0000 0 

.0000 



0.0000 

0.0000 

0.0000 0 

.0000 

B = 


0.0000 

0.0000 

0.0000 0 

.0000 


- 

0.0022 

0.0000 

-0.0388 0 

.0000 


13.5934 

0.0000 

-1.4674 0 

.0000 



0.1488 

0.0000 

-4.5577 0 

.0000 


L 

0.0000 

0.0000 

0.0000 0 

.0000 

The open- 

loop eigenvalues of the aircraft are 


The desired eigenvectors based on desired 
decoupled aircraft response are specified as 
follows, in the same order as the preceding 
eigenvalues. 5 


“v"l 0 0 

a 0 0 

q 0 0 

0 0 0 
h 0 0 
6 1 X 

p 0 0 

r X 1 
_<)>_ (1^0 

Dutch 
rol 1 


0 0X1X0 
X 1 X 0 0 0 

1 X 0 X X 0 
X X 0 X X 0 
0 0 1X10 
0 0 0 0 0 0 
0 0 0 0 0 1 
0 0 0 0 0 0 



Short Phugoid 
period 


0 

0 

0 

0 

0 

0 

X 

X 

1 


(x -arbitrary) 


The feedback gain, K, based on achievable eigen- 
vectors and closed-loop eigenvalues, is given by 


V 

a 

q 

-0.0000 

0.0000 

0.0000 

0.0038 

0.0192 

-0.0442 

0.0000 

0.0000 

0.0000 

_-0.0204 

4.1420 

-10.9580 
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h 

e 

0.0000 

0.0000 

2.7167 

0.3850 

0.0002 

0.0000 

0.0000 

0.0000 

-2.0771 

-2.6343 

-0.0010 

0.0000 

P 

r 


-0.1422 

-0.0147 

-0.0046 

0.0000 

0.0000 

0.0000 

-0.0483 

0.7094 

-0.0329 

0.0000 

0.0000 

0.0000 


-1.6433 

-1.6433 

-0.4088 

-0.4088 

-2.1413 

-0.0054 

-0.0054 

0.0000 

-0.0209 


+ 

+ 


+ 


+ 

+ 


1 . 7278i \ 
1. 7278i J 
3.2303i\ 
3.2 30 3i J 
O.OOOOi 
0.0394i\ 
0.0394i J 
O.OOOOi 
O.OOOOi 


Dutch roll mode 

Short period mode 
Rol 1 subsidence mode 
Phugoid mode 

Altitude mode 
Spiral mode 


The desired closed-loop eigenvalues are selected 
as follows 


- 2.0000 

- 2.0000 

-1.0000 

-1.0000 

-0.5000 

-0.0500 

-0.0500 

-4.0000 

- 0.0020 


+ 

+ 


+ 

+ 


+ 

+ 


4.0000i \ 
4. 0000 i J 
3.0000i1 
3.0000i J 
O.OOOOi 
0.5000i\ 
0.0500i / 
O.OOOOi 
O.OOOOi 


Dutch roll mode 

Short period mode 
Altitude mode 
Phugoid mode 

Rol 1 subsidence mode 
Spiral mode 


Two sets of trajectories are shown giving the 
results of the application of developed control 
laws to the aircraft. In the first instance, it 
is necessary to maintain the aircraft at trim con 
ditions of the specified flight conditions. No 
feedforward is required in this case. Figures 1 
to 4 show the state and control trajectories of 
the aircraft. The lateral states and the aileron 
and rudder inputs show no change from their trim 
val ues . 


In the second case, it is necessary to accel- 
erate the aircraft from Mach 0.8 to Mach 0.9, 
with a Mach rate of 0.01 Mach/sec, while main- 
taining aircraft altitude at 20,000 ft. Once the 
desired velocity of Mach 0.9 is reached, it is 
held constant. Figure 5 shows the variation of 
velocity of the aircraft when feedforward is also 
used. It shows an excellent tracking of the velo 
city. Figure 6 shows the corresponding variation 
in the altitude of the aircraft. 
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Concluding Remarks 


This paper has presented a synthesis technique 
that could be extended for control of an aircraft 
undergoing a specified maneuver. Eigenvalue/ 
eigenvector assignment procedure provides the 
necessary decoupling in aircraft handling quali- 
ties, while ensuring the location of eigenvalues to 
meet the specifications. The command generator 
tracking ensures the tracking of the controlled 
variables of the aircraft, as dictated by the 
requirements of a particular maneuver. The 
problem of developing suitable models to match the 
state trajectories of a given maneuver is being 
investigated further. In some cases this would 
involve gain scheduling and update of aircraft 
parameters, while ensuring that the constraints on 
the control surfaces are not violated. 
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Fig. 1 Variation of velocity with tine. 
Mach - 0.8 ; desired altitude = 20,000 ft; 
aircraft level. 
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Fig. 2 Variation of altitude with time. 
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Fig. 5 Variation of velocity for level accel- 
eration from Mach 0.5 to Mach 0.8. 
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Fig. 3 Elevator input to hold aircraft steady. 
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Fig. 4 Thrust input for holding aircraft steady. 
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Fig. 6 Variation in altitude for level 
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